On the first-order phase transitions in a bistable large-spin systems (Mn^Ac) 
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We have shown that recent reportEJ concerning the first-order phase transitions in the systems 
described by the Hamiltonian Ti. = —DS\ — H X S X is inaccurate. A kinetic numerical method for 
making calculations of the transition rate in a bistable system as a function of temperature has been 
00 ' evolved and a new critical value of the transition field (h x = 0.891) has been found. A physical 

t-H , mechanism caused transition from quantum to classical behavior has been discussed. 

PACS numbers: 75.45. +j, 75.50 Tt 
The thermally behavior of the escape rate is the main indicator for the search of a quantum phenomena in 
the large-spin mesoscopic systems. The system that attracted the most recent attention is the n Afni2 acetate 
(Mn 12 12 (CH 3 COO) 16 (H 2 0) A ) ■ 2CH 3 COOH ■ ^H 2 Ol, This crystal has been synthesized by LiM and its mag- 
netic properties was investigated at first by Sessoli et aB. The molecules of Mn\ 2 Ac effectively behave as magnetic 
clusters of spin S = 10 and are characterized by very strong uniaxial anisotropy Ti = —DS 2 Z with D ~ J175fc_B- 
The first observations of the tunneling of a magnetization in Mn\ 2 Ac were carried out by Novak and SessoM. They 
found that the relaxation rate of Mn\ 2 Ac followed the Arrhenius law, V ~ exp(— U e f / /T), where U e f / is height of the 
energy barrier between the two metastable states ±S, but for the certain values of the external longitudinal field the 
peaks raised. Such a peaks were. interpreted as the result of the resonant thermally assisted tunneling between the 
levels near the top of the barrierQ. Later, this phenomena attracted great experimental and theoretical interest. Th« 
successful dynamical theory of thermally activated tunneling in Mn\ 2 Ac was proposed by Garanin and ChudnovskyU. 

However, the crossover from thermal to quantum regime for the escape rate has not yet been satisfactorily explained. 
Thermal activation and tunneling are usually considered as competing processes; as a result, one might be lad to think 
i^i - that discussed mechanisms works in different temperature regions: tunneling would be observable only at extremely 
low temperatures, whereas at higher temperatures tunneling would be suppressed by the thermal activation. Using 
t-H , this approach one can divide the temperature dependence of the escape rate onto two regions - high-temperature 
region, where escape rate follows the Arrhenius law, T ~ exp(—U e ff/T) , and low-temperature region, where the 
transitions is purely quantum, T ~ exp(— B), with B independent on T. Formally, one caa,to equate the right-hand 
parts of this equations and found in such manner transition temperature T c ~ U e ff/Bu. In some sense, such a 

[ temperature may be called as a temperature of the "phase transition" from classical to quantum,behavior. Whereas 

the two -regimes smoothly join, one can call this situation the second-order phase transition!]. Chudnovsky and 
Garanina recently made an attempt to describe quantum-classical transitions in large-spin system more detailed. 
O^ ■ They consider a spin Hamiltonian 

H = -DSl - H X S X (1) 

s ■ 

with S » 1. Such a model is a good approximation for the molecular magnets Mn\ 2 Ac. The states ±S are 
separated by the energy barrier U = DS%. Overcoming this barrier is possible in the following ways: Jay absorption 
of an external field's energy, or due to thermal activation, or via tunneling. Further, authors of Ref.lH reduced the 
2 ■ initial spin problem to the particle problem using the mapping method]. The equivalent particle Hamiltonian will 
be 



H = ~+U(x) (2) 



with 



2m 



U(x) = (S + 1/2) 2 D(hl sinh 2 x - 2h x cosh a;), (3) 
and m — , h x = ^s+i) ' 
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The function U(x) in (3) described typical double- well potential (see Fig. 1) for all h x : < h x < 1. Only in 
boundary points h x — 1 and h x = this potential change the form. At h x — 1 the barrier completely disappeared, 
and when h x — > CL the width of barrier tends to infinity. 

But authors ofEJ used for the further analysis instead of U (x) its expansion with fourth-order accuracy 

V{x)=-2(s+^\ Dh x +(s + ^J D(hl-h x )x 2 + 

S+1 2}' D (l h l--^ h *)x 4 + 0(x 6 ). (4) 

For such a potential, with lowering of the transverse field below the critical value of h x = 1/4 the fourth-order term 
change the sign and became negative, and the potential "overturned". The barrier between opposite metastable 
states completely disappeared in this case, and what is more, metastable states disappeared too (see Fig. 2). It is 
difficult to believe that such situation corresponds to any physics. In reality, it means that approximation U (x) by 
V(x) is invalid. Authors of Rem alleged that field h x = 1/4 corresponds to the first-order phase transition between 
quantum and classical types of the escape (see also Ref.S). But it is easy to see, that truncation of the high-order 
terms in (4) leads to losing of the physical information. With an increasing of the accuracy of the expansion, the 
magnitude of a "critical field" rapidly decreased: if one take into account six-order term, then critical value will be 
h x = 1/24 etc. In the limit we have the exact potential (3) and the critical field h x = 0. 

Thus, it is clear that potential (4) is inadequate to the discussed problem and complete analysis of the problem for 
the initial potential (3) is necessary. For the solving of the described problem we propose the next numerical scheme. 

Let us consider a particle in potential (3). We assume that particle is in thermal equilibrium with the bath. We 
consider an ensemble containing N such particles. When the energy of the particle E less then Uq it will be localized 
in the well into one of the stationary states E n , otherwise, it will be unbound. Physically, first case corresponds 
to "pinned" magnetization, and, at last case the escape occurred and magnetization oscillated between z and z 
directions with a frequency 

7T 

^ext = 



rn f x - 
2 Jxi 



%2(E) dx 



( E ) s/E-U(x) 



where Xx^ (E) are turning points for the particle oscillating inside the potential U(x). For the values E < Uq the 
escape can took place only via tunneling. We will try to determine the fraction of the particles from the ensemble 
N which overcoming the barrier for the given temperature T due to both the thermal activation and the tunneling. 
The number of particles overcoming the barrier due the thermal activation will be 



Ntherm — N 



f°° 2 E 

/ —{k B T)- 1 - 6 VEe X p{— r -=)dE, 



where Uq — U max — U m i n , U max = —2S 2 Dh x , I7 m i n = —S 2 D(1 + h x ) . In numerical calculations, for the upper bound 
we of course make the substitution oo — > U r where U r was adjusted so as to lead to the neglect of number of particles 
outside the interval for each given value of the temperature. Residuary number of particles N we u will distribute 
among the stationary levels in according with 

exp(-Ej/k B T) 

m = N we ll (5) 

where Z — ^ exp(— Ei/k B T) is the partition function. For each Ei we calculate the probability of the tunneling 

i 

using the WKB approximation for the lowest levels 



i , i r 4 * ,„ i 



2 7->/'l,2 , • i 2 



fii = -exp(-— / {{S + -yD(h x smh ' x - 2h x cosh x) - EAdx, (6) 
7r 2D J 2 



and for the levels sear the top of barrier we use the approximation of U(x) by the parabolic potential; in this case 
probability will becl 
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Mi = t=-. (7) 

1 + eyL-pt-ZirEt/^/W) 

Ai we found numerically from the equation U(x) — Ei. The number of particles overcoming the barrier via tunneling 
will be Ntun = ^2nifli and Ntot = Ntun + Ntherm will be total number of the particles which overcome the barrier, 

i 

and in conclusion, the effective barrier transparency P one can to determine as P — N tot /N. 

For the calculation of the spectrum of the particle we use approximation U{x) by well-known potential 

W{x) = C{e- 2ax - 2e' ax ) (8) 

and assume the barrier width for the U{x\ tends to infinity (see Fig. 3). In such unperturbed case the spectrum of 
the particle inside one of the wells will bdlj 



Ei = -A I 1 - ajDjMi + I)) , (9) 

where prefactor A w I.IUq and a changes from 1.3 for h x — 0.1 to 2.4 for h x = 0.8. In real potential U(x) another 
well of course will perturb the spectrum. But the level's splitting will be essential only for levels near the top of 
barrier, which has very weak population in discussed temperature region (up to IK). Thus, for the preliminary 
estimation this approximation seems to be admissible. Let us note that in any case there exist minimal depth of the 
well E/omin- If Uq < Uq m fn-| discrete spectrum inside the well will be lacking. For the potential (8), one can found 
C^Omin from the conditionlMJ A = D/4, that corresponds to critical value of the transverse field h xcr — 0.891. It is 
clear, that for h x > h xcr the escape may occur only due to thermal activation. Therein, one can say about transition 
from quantum to classical behavior. Such a transition will be abrupt and may be called "first-order" phase transition, 
but any analogies with the Landau theory of phase transitions not stand up to criticism. 

In conclusion, let us discussed some results of the numerical calculations. From the data in Table 1 one can see 
that the spectrum of the particle in the potential (3) is characterized by the small number of levels. In discussed 
temperature region for small h x all the particles condensed, as a rule, at the ground-state level, where the tunneling 
probability is negligible. The population of the higher levels tends to zero, so, the effective transparence of the barrier 
tends in this case to zero too. With increasing of h x the population of higher levels where tunneling transparence is 
considerable, increased too, and effective transparence amount to appreciable values. For the values h x > 0.891 the 
overcoming the barrier can be possible only due to thermal activation. In Fig. 4 and 5 we plotted dependences of 
the effective barrier transparence on T for different values of the transverse field h x . The difference in the behavior 
of the curves for h x < h xcr and h x > h xcr is evident. 

Table 1. 

Some characteristics of the escape process at T — IK for the different values of the transversal field h x (N = 100). 



h x 


0.1 


0.3 


0.5 


0.6 


0.7 


0.891 


Number of levels 


7 


5 


3 


2 


1 





Population of the lowest level 


100 


100 


99.993 


99.85 


99.61 




Population of the highest level 





1.615 • lO" 11 


1.871 • 10 - 


4 0.147 


99.61 




Tunneling probability for the lowest level 








1.344-10- 


5 3.426 • 10- 


3 0.334 




Tunneling probability for the highest level 


0.261 


0.318 


0.325 


0.329 


0.334 




Probability of the thermal activation 


2.23-10" 20 


3.15-10- 11 


3.65-10- 5 


2.80-1Q- 3 


0.39 


0.62 



Figure captions 

Fig. 1. Reduced effective potential U(x) for the spin system Mn 12 Ac. h x = 1 (1), 0.75 (2), 0.5 (3), 0.25 (4), 0.15 
(5). 

Fig. 2. "Truncated" potential V{x). h x = 1 (1), 0.75 (2), 0.5 (3), 0.25 (4), 0.15 (5). 
Fig. 3. Potential U(x) and approximating potential W(x) for h x = 0.1. 

Fig. 4. Dependence of the effective barrier penetrability on T for h x — 0.7 (mixed thermal activated and tunneling 
escape). 

Fig. 5. Dependence of the effective barrier penetrability on T for h x — 0.891 (pure thermal activated escape). 
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